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Abstract 

The 4-loop sunrise graph with two massless lines, two lines of equal mass M and a line of mass 
m, for external invariant timelike and equal to is considered. We write differential equations 
in a; = m/M for the Master Integrals of the problem, which we Laurent-expand in the regularizing 
continuous dimension d around d = 4, and then solve exactly in x up to order {d — 4)'^ included; the 
result is expressed in terms of Harmonic PolyLogarithms of argument x and maximum weight 7. As a 
by product, we obtain the x — 1 value, expected to be relevant in QED 4-loop static quantities like the 
electron {g ~ 2). The analytic results were checked by an independent precise numerical calculation. 
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1 Introduction 

This paper is devoted to the analytic evaluation of the Master Integrals (Mi's) associated to the 4-loop 
sunrise graph with two massless lines, two massive lines of equal mass M, another massive line of mass m, 
with m 7^ M, and the external invariant timclikc and equal to m^, as depicted in Fig. ^ 




Figure 1: The considered 4-loop sunrise graph. 

We will follow for the analytic integration the differential equation method already proposed in [1], 
further developped in [2] (and then used in [3], [4] and in many subsequent applications which would be 
too long to report here), as well as the finite difference method [5], [6] for an independent numerical check. 

The differential equation method was already followed in two similar two- and three-loop calcula- 
tions [7], [8]; the fact that its use could be extended without major changes to the present four- loop 
calculation witnesses for its generality and power. Among the advantages of the method, it allows a rather 
clear separation of the merely algebraic part of the work (which is, not surprisingly, always very heavy in 
this kind of multiloop calculations, and can be most conveniently processed by a computer algebra program, 
in our case FORM [9]), from the really analytic issues of the problem, which can then be better investigated 
without the disturbance of the algebraic complexity. In our case, indeed, the heart of the analytic calcula- 
tion was the study of a homogeneous fourth order differential equation, whose solutions turned out to be, 
in a remarkably simple way, either a rational fraction or repeated quadratures of rational fractions. The 
required four-loop integral could then be obtained almost mechanically by repeated quadratures in terms 
of Harmonic PolyLogarithms [10]. 

Several other different approaches to the analytic evaluation of multiloop integrals are available in the 
literature, such as the powerful asymptotic expansion method (a fairly complete account can be found in 
the recent book [11]), and it would be interesting to compare the advantages and drawbacks of the various 
methods for the exacting, four-loop integration which we consider; but as the results of the present paper 
are new, a meaningful comparison cannot yet be carried out. 

Following, as already said above, the approach already used in [7] and [8], we identify the Mi's of 
the current problem within the continuous d-dimensional regularization, write the system of differential 
equations in a; = m/M satisfied by the Mi's, convert it into a higher order differential equation for a single 
MI, Laurent-expand in [d — 4) around d = 4, solve the associated homogeneous equation at d = 4 (as in 
previous cases, the solutions of the homogeneous equation are surprisingly simple) and then use recursively 
Euler's method of the variation of the constants for obtaining the coefficient of the [d — 4) expansion in 
closed analytic form. The result involves Harmonic PolyLogarithms (HPL's) [10] of argument x and weight 
increasing with the order in (d — 4). We push the analytic calculation, which works up to virtually any 
order in (d — 4), up to {d — 4)"^ included, involving HPL's of weight up to w = 7 included. The integration 
constants are fixed at a: = ; as a by product we obtain the values at x = 1, which are relevant in the 
evaluation of 4-loop static quantities such as the electron [g — 2) in QED. The result was checked and 
confirmed by an independent numerical calculation performed with the method of [5] , [6] . 

The plan of the paper is the following. In Section |21 we define the Master Integrals and write the 
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differential equations; in Section |31 we study the x behaviour; in Section^ we work out the Laurent- 
expansion in (d — 4) and discuss the associated homogeneous equation; in Section [S] we write down the 
solutions by Euler's method of the variation of the constants and give the x = 1 values; Section deals 
with the independent numerical calculation by which we check the results at a; = 1 of the previous Section, 
in Sectional we carry out the evaluation of the a; — > values by direct integration in the parametric space. 



2 The Master Integrals and the Differential Equations 

We find that the problem has 5 Mi's, which we choose to be 

1 f d'^ki d'^k2 d'^k^ d'^k^ N, 



(27^)4(^-2) J klkl{kl + A.P){kl+A.P)[{P -ki-k2~k:i~-kiY + m^] ' 

(2.1) 

where the 5 numerators Ni are (M^, fci • fc3,p • ^3, fci • fc2,p • ^2)- In terms of the dimensionless variable 
X = m/M and putting P — Mp one can introduce 5 dimensionless functions ^i{d, x) through 

F,{d,M^,m'^,P^ = -m^) = Af^'^"* C^{d) $,(d,a;) , (2.2) 

4 — d 

where C{d) ~ (47r)~2-r(3 — d/2) is an overall loop normalization factor, with the limiting value C(4) = 1 at 
d — A. Some of the formulae which will follow (in particular the differential equations) are slightly simpler 
when written in terms of x"^ rather than x; but as x is the most convenient variable for expressing the final 
analytic results, we stick to x from the very beginning. 

As in [7] [8], the derivation of the system of differential equations is straightforward; the derivatives of 
the Mi's, i.e. of the 5 functions $i((i, x), with respect to x are easily carried out in their representation 
as loop-integrals Eg. 1)2. II) : when the result is in turn expressed in terms of the same Mi's, one obtains the 
following linear system of first order differential equations in x 

d^i(d,x) f(3rf-7) 3(rf-2) 3 (d - 2) 1 ^ , , , [3(d-2) 3 (d - 2) 

- ' ''9i[d,x) ' 



d^2{d,x) 




dx I X 2{l-x) 2{l + x) \ ^' ' 1 X 2{l~x) 

3^2{d, x) - 3$3(ci, x) + $4(d, x) - 3$5(d, x)] , (2.3) 



dx 

d<^3{d,x) _ j 3 (d - 2) 3 (d - 2) 
di ^ "I2 (1 -x) ~ 2~(1 

f3(d-2) 9(d-2) 9(d-2 



$2(rf, x) - 2$3(d, x) + Md, x) - 2$5(d, x) , (2.4) 



1 X 2 {1~ x) 2 (1 +a;) 

d<^i{d,x) _ 2{d-2 



^i{d,x) - 3$2(rf, a;) - $4(d, x) 

^3{d,x)+Md,x)) , (2.5) 



dx x 
d^^{d,x) 2{d-2) 



$2(d,x) + $4(rf,x) , (2.6) 



dx 



$3(d,x) + $5(d,x) , (2.7) 



At variance with the cases discussed in [7], [8], the system is homogeneous; indeed, quite in general the 
non homogeneous terms are given by the Mi's of the "subtopologies" of the considered graph, obtained 
by shrinking to a point any of its propagator lines. When that is done for the 5-propagator "topology" of 
Fig-O; one obtains the product of 4 tadpoles; but as the considered graph has two massless propagators, at 
least one massless tadpole is always present in the product; as in the d-dimensional regularization massless 
tadpoles vanish, the product of the 4 tadpoles is always equal to zero - and therefore the differential 
equations are homogeneous. 
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By inspection, one sees that $3(d, x), $5 (d, x) appear in the r.h.s. of Ea.s (|2.3l2.7|l only in the combina- 
tion 

*3(d, x) = $3(d, x) + $5(d, x) ; (2.8) 
the other hnearly independent combination of the two Mis, say 

^5id,x) = ^3{d,x)-^5{d,x) , (2.9) 

decouples and can be expressed in terms of the other integrals by means of the trivial 1st order differential 
equation 

d'^sid, x) _ [ 3 (rf-2) 3 (d- 2) ' 
dx ~ "[2 (1 -x) ~ 2 (l + x) 

f5(d-2) 9(d-2) 9{d-2) 



$i(d, x) - 3$2(d, a;) - $4(^,2;) 

*3(d,x) (2.10) 



1 X 2(l-a;) 2(1 + 2;) 

As 'i>5{d,x) does not enter in the r.h.s. of Ea.s (|2.3l2.7|l the 4 linear equations for <i>i (d. x). ^2{d,x), ^3((i, x) 
and $4((i, x) can be written as a fourth order equation for the first Master Integral ^i{x), which will be 
called simply $((i, x) from now on, and which is therefore equal to 

._CzHd^[ d'^fci d'^k^ d'^k, d'^h (2__2. 

^ ' ~ (2^)4('i-2) J klklikl + l)(fc| + l)[(p -k,-k,-ks^ k,Y + X2] ' - ^ ^ ■ ^^-''^ 

One obtains for $((i, x) the following 4-th order differential equation 

x^(l - x^)'^^ + x^{l + 5x^ - 3(d - 4)(1 - 3x^)}^^^^ 

d2$(d,x) 



-x|l2 -I- 6x2 „ 4)(i3 32a;2) + (^^ _ 4)2^^ _^ 26x2) 
-H<{ 12 - 18x2 + [d- 4)(25 - 2x2) + 8(d - 4)2(2 + 5x2) _^ _ 4^3^-^ ^ g^2^ 



(ix2 
(i<i>(d, x) 



c?x 

Axl + 12 + 29(d - 4) + 23(d - 4)2 + 6(rf - 4)^ U(d, x) = . (2.12) 
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The expression of the other Mi's in terms of <i>(d, x) and its first 3 x-derivatives reads 



11 x'^ 2x'^ 7x'^ 

a;(l - 7x2) '^Q^2-^ ^(^11 28x^) \ x) 



(x{l-7x^) x{l + 16x^) x{n + 28x^)^( 
1 10 (d-2)3 ^ 12 (d-2)2 ~ 60 (d- 2) j 
a;2(l-2a;2) a;2(l - lOx^) ] ^^^(d, x) 



da: 



5 (d - 2)3 30 (d - 2)2 J dx2 
a;3(l-a;2) d3$(d,a;) 
~20(d-2)3 da:3 ' 

^^^'^'^^ = -{ 15 (d - 2)3 - 3F^ + 

a;(l-7x2) a;(l-8a;2) a:(l - 12^2) 1 d$(d, x) 
30 (d- 2)3 ~ 12 (d- 2)2 ^ 20 (d - 2) J " 



(2.13) 



x2(l-2x2) j;2(2-5x2)\ d2$(d,x) 



da; 



15 (d ~ 2)3 30 (d - 2)2 J dx2 
a;3(l-a:2) d3$(d,x) 
~60 (d- 2)3 dx3 ' 

^ .J ^ f 1 1 4x2 43^2 

'^^('^'"^ = -l2-6(d^ + ^Td3^-(d^+15(d^!'^('^'"^ 



(2.14) 



2x x(14-3Gx2) x(ll + 43x2)1 d$(d,x) 

5 (d-2)3 " 5 (d-2)2 15 (d - 2) J 

4x2(1-2x2) 3x2(1- 5x2) \d2$(d,x) 



dx 



5 (d - 2)3 10 (d - 2)2 J dx2 

x3(l -x2) d3$(d,x) 



(2.15) 



5 (d - 2)3 dx3 

3 The X ^ behaviour of ^{d, x) 

By inspection, one finds that the most general solution of Ea. ()2.12|) can be expanded for x ^ in the form 



$(d,x)=^x"' (^AW(d)x2") , (3.1) 
where the values of the 4 exponent ai are 



\n=0 



ai = , 

a2 - (d - 2) , 

as = -id -2), 

ai = {M-7) ; (3.2) 

the A'^\d) are the 4 integration constants, and all the other coefficients An\d) for n > are determined 
by the differential equation Ea. 1)2. 12(1 once the integration constants are fixed. 

It is interesting to recall the leading exponents of the x — > expansions of the solutions of the corre- 
sponding equations encountered in [7] and [8]. Calling a] the exponents for the 2-loop graph of [7] and 
al those of the 3-loop graph of [8], one finds that there are always 4 exponents. The explicit values in the 
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case of [7] are 



v(2) 







(2) 
"3 



= -id-2), 

af^ = (d-3) (3.3) 



and in the case of [8] 



= 0, 



J3) 



(d-2) , 
-{d-2) , 

(2d -5), (3.4) 



where the exponents marked by a bullet (•) correspond to the behaviours forced by the inhomogcneous 
terms. The similarity between the 3 sets of behaviours is impressive: the first 3 exponents are identical, 
the fourth differ in steps of {d — 2) for each additional loop. 

In more details, 2 of the exponents of Ea.s (|3.3|l correspond to the 2 independent solutions of the 
associated homogeneous differential equation, which is of 2nd order, while the other 2 are forced by the 2 
independent behaviours for a; — > developed by the inhomogcneous terms (both products of 2 tadpoles, the 
first product of two tadpoles of mass M, the second of a tadpole of mass M and a tadpole of mass m = Mx). 
In the a; — > expansion of the most general solution, one is therefore left with 2 undetermined integration 
constants, corresponding to the two homogeneous solutions, while all the terms with the behaviours of the 
inhomogcneous terms are fully determined by the inhomogeneous equation itself. Similarly, in Ea.s H3.4(l 
3 exponents correspond to the 3 solutions of the 3rd order homogeneous equation, with the 4th exponent 
forced by the inhomogeneous term (product of 3 tadpoles, of masses M,M and m = Mx), so that there 
are in principle 3 undetermined integration constants, the term corresponding to the remaining behaviour 
being fixed by the equation. In the present case, finally, the equation is homogeneous (as already observed, 
there is no inhomogeneous term, as all the possible products of 4 tadpoles involve at least a vanishing 
zero-mass tadpole) and the 4 exponents of Ea.s (|3.2|l correspond to the behaviour of the 4 homogeneous 
solutions, so that one is left in principle with 4 undetermined integration constants. 

A qualitative inspection of the integrals which one tries to evaluate by means of the differential equations 
fEa. (|2.11|) of the present paper, Eq.(1.3) of [7] and Eq.(2.3) of [8]) shows that they are all finite (just finite, 
not analytic!) for x — > 0+ and {d — 2) > 0; that is sufficient to rule out from their expression as solutions 
of the differential equation the terms with the behaviour of the third and the fourth exponent (which is 
negative when d is just above 2). 

In the case of [7], that fixes completely the solution. In the case of [8], one integration constant is left 
undetermined; to fix it, one has to provide some independent information, such as the value of the required 
Fcynman integral at a; = (which corresponds to a simpler vacuum amplitude); that value can be provided 
by an explicit "conventional" calculation, say in parameter space, which is in any case much easier than a 
calculation for non-zero values of the variable x. (A closer analysis carried out in [8] shows however that the 
regularity of the integral at x = 1 is sufficient to fully determine the solution, so that the actual knowledge 
of the X = value can be used as an independent check). 

In the current case, as the equation for ^{d,x) is homogeneous, the only information is that (d) 
and A'^'' (d) are both equal to zero due to the finiteness for a; ^ 0+ ; by substituting the ansatz Eq. ()3.1|l in 
Eq.(E13l and dropping a';^^ {d),A'i^^ (d), one finds for $(d, x) Ea. lfTTTl) the a: expansion 

- 2(2d-5)(3d-8) , ^,^4. 
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The expansion depends on the two as yet unspecified integrations constants {d),AQ {d) - which are 
to be provided by an independent, expUcit calculation. That is done in Section [7| see Ea. (|7.14|l . by direct 
integration in the parametric space. 

Let us note here, for completeness, that in the present case the knowledge of the regularity of the 
solution at X = I does not provide any additional information. 



4 The expansion in ((i— 4) and the homogeneous equation at d = 4. 

The Laurent's expansion in (d — 4) of ^{d, x) Ea. (|2.11|l is 



(4.1) 



n=— 4 



as it is known on general grounds that it develops at most a fourth order pole in (c? — 4). By substituting 
in Ea. H2.12|l one obtains a system of inhomogeneous, chained equations for the coefficients $^"^(a;) of the 
expansion in (c? — 4) ; the generic equation reads 



dec* dx'^ dx^ dx 



$(")(a;) =i^(")(a^) , (4-2) 



with 



i^(")(a;) = <j 24a; + (3 + 24x2)— 



d 



(x) 



92x + (16 + 40x2) — - {x + 26x^ 
dx 



+ (ll6x+(25-2x2)A_ (13:^ + 32x3)-^ - (3x2 -9x4)-^U("-i)(x) ^ (43) 
^ ' dx dx^ dx-^ J 

which shows that the equation at a given order n for $'^"^(x) involves in the inhomogeneous term the 
coefficients $'^'^)(.x) (and their derivatives) with k < n - hence the "chained equations" expression used 
above (obviously $'^'^)(x) = when k < —4). Such a structure calls for an algorithm of solution bottom-up, 
i.e. starting from the lowest value of n (which is n = —4) and proceeding recursively to the next n + 1 
value up to the required order. 

The Ea.s (|4.2(l have all the same associated homogeneous equation, independent of n, 



xHl - + + 5x2)-^ - 6x(2 + x2)-^ + 6(2 - 3x2)-^ + 48x 

dx* dx'^ dx^ dx 



0(x) = 



(4.4) 



once the solutions of Ea. H4.4|l are known, all the Ea.s (|4.2| l can be solved by the method of the variation of 
the constants of Euler. 

To our (pleasant) surprise, the solutions of Ea. l4.4|l are almost elementary. By trial and error, a first 
solution is found to be 

0i(x)=x2. (4.5) 



We then substitute (^(x) 
derivative of ^(x) 



i(x)^(x) in Ea. H4.4|l . obtaining the following 3rd order equation for the 

Cix) = , (4.6) 



x3(l - x2)-^ + 3x2(3 -x^)^+ 6x(l + 2x2)-^ - 6(5 + 2x2) 
dx'' dx^ dx 



and find that it admits the solution 



1 



C^(x) = ^(l-x2+x4) 



(4.7) 



6 



Substituting £,'{x) — £,'2{x)x{^) in Ea. (|4.6|) we obtain the following 2nd order equation for 



W2 

x^(l - x^)(l - + x")— + 6x^(2 - x^)— - 6(2 - 2x'^ + 
dx-^ 



d_ 

' dx 



X'{x) = 



which admits as solution 



5~~2x^ + 5a 



x-^ (1 — + x^y 



(4.8) 
(4.9) 



Finally, substituting x'(a;) — X3{x)t{x) in Ea. (|4.8|) . we obtain the equation 

x{l - a;2)5(l - x' + - 2x' + 5x^)4 

ax 

-2(1 - 2;2)4(15 - \2x^ + llx* + 30a;^ - 24x* + 2Qx^^) 
which has the solution 

, x^ 1 — x^ + x"^ 



/(x) = , 



(4.10) 



(4.11) 



(1 - x2)5 5 - 2a;2 + Sx^ 

By repeated quadratures in x and multiplications by the previous solutions we obtain the explicit 
analytic expressions of the 4 solutions of Ea. (|4.4(l : the nasty denominators {1 — x'^ + x^) and (5 — 2x^ + 5x'*) 
disappear in the final results, while the repeated integrations of the terms with denominators x, (l + x) and 
(1 — x) generate, almost by definition, Harmonic PolyLogarithms [10] of argument x and weight increasing 
up to 3. Explicitly, one finds 



(t>2ix) 
04 (x) 



1 



(1 -x") -H(0,x)x2 



(5 + 18x2 + 14x^ + 5x**) 1 



(12 + x^ - l2x^)H{0; x) + 12 x^H{0, 0; x) , 



8 x2 2 

(1 + x2)(15 + 182x2 + ISx'^) ^ 3(1 - x2)2(5 + x2)(l + Sx^) 
65536 X 

9(1 -X4) r 



8192 



131072 x2 
9x2 , 



H{-l;x)+H{l;x) 



4096 



7J(0,0, -l;x) + i7(0,0,l;x) 



(4.12) 
(4.13) 

(4.14) 



W{x) 



(1-X2) 



2^3 



(4.15) 



H(0,-l;x) + i7(0,l;x) 
The corresponding Wronskian has the remarkably simple expression 

4>l{x) (t>2{x) (t>z{x) (t>A{x) 
C^'^ix) 0'2(x) 0(,(x) 0Ux) 

<(x) 0^'(x) 0(((x) <jy'i{x) 

0;"(X) 0-(x) 4>'i'{x) 4>'i'{x) 

in agreement (of course) with the coefficients of the 4*'' and 3'"'' x-derivative of (/)(x) in Ea. H4.4|l . 



5 The Solution of the differential equations for the coefficients of 
the expansion in [d — A). 

With the results established in the previous Section one can use Euler's method of the variation of the 
constants for solving Ea.s H4.2|l recursively in n, starting from n — —4. We write Euler's formula as 



$(")(a;) = (t>,{x) 



i=l 



(«) 



dx' 



W{x 



-Nh{x')K'^^\x') 



(5.1) 
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where the (jjii^) are the solutions of the homogeneous equation given in Ea.s ()4.5l4.14(l . the $1"' are the as 
yet undetermined integration constants, the Wronskian W{x) can be read from Ea. H4.15|) . the AIi{x) are 
the minors of the (j)^ (x) in the determinant Ea. H4.15|l . and the K'^^^x) are the inhomogeneous terms of 
Eq. (|4.3|) . The constants are then fixed by comparing the expansion in x for x — > of Eq. (|5.1|) with 
the expansion in {d — 4) for d — > 4 of Eq. H3.5|l . Exphcitly, we find for instance 

1 o 



= --X-, 

^ ^ 384 256 192 48 ^ ' ' ' 



(5.2) 
(5.3) 



$(-2)(x) = 
$(-i)(x) = 



41 



143 



4608 3072 



37 
2304" 



9216 



3 1 

64 ~ 96" 



x' ] x^H{0; x) - j-x'^H{0, 0; x) , (5.4) 



48 



805 



205 



55296 4096 



1 



743 4 599 

rX + , . . X 



27648 



,64 96 
3173 1 



221184 

1 



43 



11 



768 384 



— - — x^ x^i?(0,0;x) - —x''H{0,0,0;x) , 



48 



663552 3840 
5 1 



4 2 
-TT X 



15649 



49152 96^^^^ 331776^ 



2304 



5308416' 96^^ ^ 



(5.5) 



1536 48 



C(3)x^ 



1 

768" 



121 
2304" 



599 
55296" 



x" ]H{0;x)+ — + 



1 



256 



x2 + 



1 

96' 



1536 



i7(l,0;x) -i7(-l,0;x) 



43 
768^ 



5 1 



96 1536x2 
11 



384 



+ 2304^ + 



-iJ(0,0,-l,0;x) 



32 
_3_ _ J_ 

64 " 96" 



16 



iJ(0,l,0;x) -i7(0, -l,0;x) 

x2 ) x^H{0, 0, 0; x) - -^x^H{0, 0, 0, 0; x) 
48 



i?(0,0,l,0;x) 



(5.6) 



The full results become quickly too lengthy to be reported explicitly here, so we give only the values of the 
integration constants up to order 3 included. We find $4*^' = 0, for (all) — 4 < fc < 3; the other constants 
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are 



$W = 24978775 ^241^ ^18^^, _^ 

^ 603979776 786432^^ ^ 5242880 64^^ ' ^ ' 

= _i90M2^^^^(3)_^50^,4^ (5.3) 

^ 905969664 393216^ ' 7864320 ^ 

= 264575 _^^(3) ^^.^ (5.9) 

^ 452984832 196608^^ ' 3932160 ^ ^ 

= 1791217627 _ i76887 37601 ,4 ^ ^^(5) ^ 1^.(3) _ ^ (5.IO) 

^ 7247757312 3145728^^ ^ 62914560 262144^^ ' 32^ ^ ' 32256 ^ ' 



^(.) = 0-4733523 _ ^9im 79009 _ ^ 

^ 10871635968 4718592^^ ^ 18874368 131072^^ ' ' ^ ' 

^ _ 191322157 J704^ _ ^78^ _ 

^ 5435817984 2359296^^ 9437184 65536^^ ^ ' 

^(3) _ _ 40220928899 6041939 775961 4 3855 AC5^_ 12543 a.gx 

^ ~ "21743271936"'' 9437184^^ -* " 150994944 1048576^^ ' ~ 131072^ *■ ' 
54353 fi 1 4 605,,„, 

+ 44040192'^ + SSi^^)'^ - 5T2«^) ^ ^'-''^ 

^(3) 336259373 13GG211 , ^ 7536877 , 1939871 , , 15281 .o.ox 

*9 = 1 (, 3 7r*H C(5) H C (3) 

^ 16307453952 14155776^^ ^ 1132462080 1572864^^ ' 196608^ ^ ^ 

TT^ : (5.14) 

66060288 ' ^ ' 

(3) ^ 3741658007 381061 275627 41641 19 

^ 16307453952 7077888^^ 566231040 786432^^ 32768^ ^ 

247 

(5.15) 



33030144 



Figure 2: The on-shell 4-loop sunrise graph. 
As a byproduct, using the results of [12], we obtain the values at a; = 1, i.e. the on mass-shell values 
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at —p^ = 



M2 of Eq.lEU, depicted in Fig.2, 



= 1) 



1 



7 



17 



835 



7379 



1 



64 {d - Ay 



256 (d - 4)^^ 
53 



, , 6766055 
^^'-'^^ - 42467328 + 
1449210865 



768 (d - 4)2 

9 . 7 A 



73728 (d- 4) 1769472 1920 



+ (d-4)^ 



+ (d-4)= 



1019215872 
203 



1152 

1913^ 

4608' 



7680 
53 



128 



7r^ln2 



7 
256 
1855 



C(3) 



29 
256 



C(5) 



17 



2304 



23040 



49 
1024 



^'C(3) 



1024^^ ^ 256 ^ 32 256 256 



245 



rC(3) 



661 



1024' ' ' 322560 
182188906799 734695 



24461180928 
8707 



147456 
21 



265, . 

ln*2 

768 



331776 
265 



32 



■a4 



1913 
512 
119 



7r2ln2 



1024 



7r2ln42 



21 



-— 7r2C(3)ln^2 
512 ' 

-^A(3) + 
12288 ' 



1024 
2. 123 



3072 
661 



^'C(3) 



66955 
9216 
493 



3072 



1280 



184320 
9917 



4096 
45 



128 



^C'(3)ln2 + — C(3)ln42- 
1024^ ^ ' 128^^ ' 

81 2, r. 45 81 2 45 , 2„ 

^32'^'*'^ + Y^'^4C(3) - ^"s'^ ^ Te"^ 

45 , „ 45, , „ 135 45, 45 , 
-— a6ln2 + Y^06ln2 —aj + yOy - + 



a47r - 

107 

1536 
15 

' 2048 
1395 



C(5)4 

1715 
4096 



583 
384 
441 



C(3) + 



1024 
C'(3) 



7r2C(3)ln2 



7r4ln32 



7rSln2 



1024 
15 



((5)111^2 

29335 



1792 



ln^2 



2048 



<(7) 



+o{{d-iy 



(5.16) 



At variance with [12], we have expressed the results in terms of the constants hsted in Table 1; the first 
column is the name of the constant, the second column its value as Harmonic PolyLogarithm of suitable 
argument, the third as Nielsen PolyLogarithm (when available), the last the numerical value. 



constant 


HPL 


NPl 


numerical value 


C(3) 


i/(0,0,l;l) 


^2,1(1) 


1.2020569031595942854 




i/(0, 0,0,1; 1/2) 


^3.1(1/2) 


0.51747906167389938633 


C(5) 


i/(0,0,0,0,l;l) 


^4.l(l) 


1.0369277551433699263 


as 


i7(0,0,0,0,l;l/2) 


^4,l(l/2) 


0.50840057924226870746 




i/(0, 0,0, 0,0,1; 1/2) 


^5.1(1/2) 


0.50409539780398855069 


be 


i7(0, 0,0, 0,1,1; 1/2) 


^4.2(1/2) 


0.0087230030575968884272 


C(7) 


i/(0,0,0,0,0,0,l;l) 


56,l(l) 


1.0083492773819228268 


aj 


i/(0, 0,0, 0,0, 0,1; 1/2) 


^6,i(l/2) 


0.50201456332470849457 


b7 


i/(0, 0,0, 0,0, 1,1; 1/2) 


^5,2(1/2) 


0.0041965726953603256975 


dj 


i7(0,0,0,0,l,-l,-l;l) 




0.0022500546439578516764 



Table 1: Constants up to weight 7 appearing in Ea. l|5.16|l 

Once the explicit analytic expressions of the coefficients of the Laurent-expansion of ^{d,x) in {d — 4) 
are known, one can use Ea.s (|2.13l2.15|) for obtaining the coefficients of $2(rf, x), 5'3((i, x) and ^^{d, x). The 
coefficients of the expansion of ^'5((i, x) are then recovered by integrating in a; Ea. H2.10() : the quadrature is 
trivial to carry out in terms of HPLs, and the integration constants are fixed by the condition 5*5 (d, 0) — 0, 
which follows at once from Ea.s (|2.9|I . H2.2|l . H2.1|l . 
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6 The Independent numerical calculation. 



In this Section we will calculate $(c?, a; = 1) with high numerical precision by suitably using the difference 
equation method described in [5,6] and references therein. At variance with [6], we will apply the method 
directly to the Master Integral with massless lines. This will imply some additional work for establishing 
the initial conditions of the difference equations. 
We define 



-2d 



d'^h d'^k2 d'^ks d'^ki 



+ l)»((fc2 - fci)2 + l){{k3 - k2Y + \){ki - k^Yip - kiY ' ^ 

so that /5(ci, 1) is equal to x) of Ea. (|2.11|l at x = 1 up to a known multiplicative factor 

h{d, 1) = [4r(l + e)f<^{d = 4 - 2e, x 1) . 



-1 , 



(6.1) 



(6.2) 



By combining identities obtained by integration by parts one finds that /s (d, n) satisfies the third-order 
difference equation 



32(n - l)(n - 2)(n - 3)(n - 3d + 5)/5(d, n) 



'4(n-2)(n-3) 



+2(n- 3) 



133d - 84)n + ~ Uld^ + 134d - 24 



h{d,n- 1) 
h{d,n~2) 



15^2 + (39 - 50d)?i + 27^2 + 5d - 54 

12n^-(38(i+ 24)^2 + (23^2 

+ {n-d- \){n -2d+ l)(2n - 3d)(2n -5d + 4,)h{d, n - 3) = 
We will solve this difference equation by using the Laplace ansatz 



(6.3) 



I^{d,n) = / r-'-V5id,t) dt 



(6.4) 



giving for {d, t) the fourth-order differential equation 

d^ 



4t4(8i-f l)(t-l)2^i;5(d,t) 
+^^{1- 1 



dt^ 

24(d -I- l)t^ -h (18 - 2Qd)t -Id +12 



dt^ 



V5{d,t) 



t 

+t 



(576(d - l)t^ + (-108^2 - 420d + 6A8)t^ + (46^^ -I- 38d - 144)i + Tld^ - 284d + 288 



dt^ 



V5{d,t) 



(576d - 960)t^ + (-216^2 + 360d)t^ 



+ {-l8d^ + 190^2 - 496d + 384)t + 77^^ - 533^^ + 1236d - 960) 
-f(d- 3)(2d- 5)(3d- 8)(5d- 12)w5(d,t) = . 



j^v,id,t) 



(6.5) 



We will look for the solution of Ea. H6.5|l in the form of a power series expansions; inserted in Ea. H6.4|l and 
integrated term by term it will provide very accurate values of Iz{d,n). As the convergence is faster for 
larger n, we will consider "large enough" values of the index n (see below); the repeated use "top-down" 
of Ea. (|6.3() {i.e. using it for expressing I^[d, n — 3) in terms of the /5(fc) with k = n,n — l,n — 2) will give 
the values corresponding to smaller indices, till I^id, 1) is eventually obtained. To go on with this program, 
initial conditions for v^^djt) are needed. 
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From the definition Eq. , and introducing spherical coordinates in d-dimension for the loop momen- 
tum fci, d'^ki — k'l^^dkidD.{d,ki) one has 



Md, kl) 



(fc2)''/2-i dkl 



h{d,kl) , 



dn{d, kl) 
n{d) 



h{d,i,{p-kif) 



(6.6) 



(6.7) 



where f2((i) is the d-dimensional solid angle, and l4{d, 1, {p — ki)^) is the 3-loop (off mass-shell) sunrise 
integral 



h{d,n, (p - ki)^) = n 



-3d/2 



d''k2 d'^ks d''ki 



((fc2 - ki)^ + - + l){ki - ksfip - k^Y 

By the change of variable l/(fcj + 1) = <, k'l = {\ — t)/t, one finds 

1 



h{d^n) = ^j e'\\-t)i-H~ih[d,^-^^ dt, 



from which one gets the relation between v^{d,t) and f^{d^ (1 —t)/t) 



V5{d,t) = (1 -f) 2-4-^/5 —^j 



dt . 



(6.8) 



(6.9) 



(6.10) 



From that relation we see that we can derive boundary conditions for (d, t) in the t —> I limit from the 
expansion of fsid, kf) in the fci ^ limit, which is easy to obtain. Indeed, only the first denominator of 
Eq. (|6.8(l depends on ki ; expanding for small fci and performing the angular integration one gets 



dn{d, fci) 



dn{d, fci) f 1 fc? - 2fci ■ fca (fc? - 2fci • fc2)2 

^ (fci + 1)^ 



n{d) (fc2-fci)2 + i J n{d) Vfcf + i (fcf + i) 



"fc^ + i 



+ kt - 



4 fci 



(fci + 1)2 d (fci + 1)3 



(6.11) 



The above result gives the expansion of /^{d, kf) at fc^ ~ 0: 

h{d,kl) =/f (d) + /«(d) kf + 0{kt) 



/f(d) =/4(d,l,p') 

/i'^(d) = -hid,2,p^) + ^[h{d,2,p^)-h{d,3,p^)] 



(6.12) 



Note that /5(d, fc^) is regular in the origin. 

By inspecting the differential equation (|6.5|l one finds that the behaviour at i = 1 of the 4 independent 
solution is ^ (1 — i)"S with ai = d/2 — 1, a2 = d/2, — 0, and a4 — 1; for comparison with Ea. H6.10(l 
the behaviours — 0, and = 1 arc ruled out and the expansion reads 



v,{d,t) - (1 - t)^-i (4°^(d) + vi'\d) (1 - t) + 0(1 - tf) ; 
by comparison with Ea. H6.12(l (t = 1 corresponds to fc^ = 0), one obtains 



(6.13) 



v'^\d)^ 



r(i)- 



1 



ri 



(6.14) 
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The values l4{d, n) of Ii{d, n^p^) at 



1 are therefore required 
d'^k2 d'^k-i d'^ki 



{{k^-k2? + l){ki-k^Y{p-k^Y- 



P 



-1 . 



(6.15) 



The problem of evaluating the Ii{d,n) is similar to the original problem of evaluating the l5{d,n), but in 
fact it is much simpler, as the l4,{d,n) involve one less loop and one less propagator. As above, by using 
integration- by-parts identities one finds that Ii{d,n) satisfies the third-order difference equation 

6(n - l)(n - 2)(n - ^)h{d,n) 

-{n - 2){n - 3)(10n -Id- I0)h{d, n - 1) 
+ (n - 3)(2n2 + (2d - 18)n - 7d^ + 2M - 8)h{d, n-2) 
+ {n-d- l)(n -2d+ l)(2n - 3d)hid, n - 3) = . 



We solve the difference equation by using again the Laplace ansatz 

1 

l4{d,n) = j t'^~'^V4{d,t) dt , 



where V4 {d, t) satisfies the differential equation 

^3 

2t^{U+l){t-lf-^Vi{d,t) 



(6.16) 



(6.17) 



+i^{t - l)(36f^ + (6 - 7d)t -M+ lS,)—,vAd, t) 

dt'^ 

+t{m^ - Udt^ + {-7d^ + 33d - 36)t + 13d^ - 61d + 72)-^V4{d, t) 



dt 



+ {d- S){2d - 5)(3d - 8)v4{d, t)=0. 
Following the procedure used above, we write 



h{d, n) 

U{kl) 
hid,{p-k2f) 

1)4 (d, t) 



hi 



(fej)*^^^ ^ dk2 f ,,2\ 



/ 



dn{d, k2) 
0(d) 



h{d,{p-k2f) , 
d'^kz d'^ki 



((fc3-fe)2 + l)(fc4-fc3)2(p-A;4)2 
1 - r 



r(i) 



{i-t)-^-h-^U(d 



t 



(6.18) 

(6.19) 

(6.20) 
(6.21) 
(6.22) 



At variance with the previous case, the function /4((i, fc|) is not regular for k2 0, as at fc2 = the value 
of the external momentum squared {p — fc2)^ becomes the threshold of the 2-loop sunrise graph associated 
to I^{d,p'^). But it is not difficult to evaluate analytically /3(d, q^) for generic off-shell by using Feynman 
parameters: 



h{d,q^) 



2r(5-d)r(3-f)r^(f-i) 
(d - 4)2(3 -d)r(f) 



,i^i(3-d,2-^;^;V 



where 2-F1 is the Gauss hypergeometric function. The expansion of /3(d, q^) in = —1 consists of the sum 
of two series, 



h{d,q^) = ao{d) 



1 + 0{q^ + 1) 



+ bo{d) (g2 + 1) 



2d-5 



1 + 0(g2 + 1) 



(6.23) 
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ao(d) 



bo{d) 



2r(5 - d)r (3 - I) r2 (| - i) r(2d - 5) 



(4-d)2(3-d)r 

2 



3)r(d-2) 



1 r(5 - 2d) 



Inserting Ea. (|6.23|l into Ea. H6.20() . setting q = p — k2 and performing the angular integration over k2 by 
means of the formula (see Eq.(88) of Ref. [5]) valid for fc2 ^ 



with N = —{2d — 5), as in the term {q-^ + 1 



r(|)r(f) 



2T{N)r{^{d- N)) 



of Ea. Ht).23|l . one obtains 



h{d,kl) = ao(d) 



1 + 0(A:2) 



r(|)r(f-d) 



2r(5-2d)r(i(3d-5)) 



1 + 0{kl) 



(6.24) 



(6.25) 



Using the variable 1/(A;| + 1) = < in Ea. (|6.25|l and inserting it in Ea. (|6.22j) one gets the initial condition 
for W4 [d, t) at the singular point t = 1 



l + 0(l-t) 



r(l-d) 



2r(5-2d)r(i(3d-5)) 



feo(d)(l-t)^(3'^-^) 



l + 0(l-t) 



By inspecting the equation H6.18|l one gets that the behaviour at i = of {d, t) is 



Vi{d,t^ 0) « ci'V) i^''^^ + ci\d) r 



2d+5 



-cf (d) t 



i(-3d+8) 



(6.26) 



(6.27) 



so that for d ^ A the integral (|6.17() is convergent for n > 4. 

All the quantities depending on d are then systematically expanded in d — 4, the series are truncated 
at some fixed number of terms, and the calculations with the truncated series are performed by using the 
program SYS [5]; as the first 12 terms of the series are lost in the intermediate steps of the calculations, 
in order to obtain the final results, from 1/e'* up to O(e^), 23 initial terms are needed. We solve finally 
the differential equation 1)6. 18(1 with the initial condition (|6.26() by a first expansions in series aX t = 1; 
due to the presence in Ea. (|6.18|l of a singular point at t — —1/3, to have fast convergence till i = we 
switch to subsequent series expansions at the intermediate points 1/2, 1/4, 1/8 and 0; then we calculate 
the integral H6.17|l for ti = 4, 5, 6, 7, 8 by integrating the series term by term (about 300 terms are needed 
to reach a precision of 77 digits). By applying repeatedly "top-down" the recurrence relation H6.16|l to 
liid, 8), lA{d, 7), l4,{d, 6), we obtain Ii{d, 5) and /4((i, 4) (which are cross-checked with the values obtained 
by direct integration), then Ii{d, 3), Ii{d, 2) and Ii{d, 1) 



/4(d,l) =r(i + e) 



0.3333333333333333333333333333333e- 



- 1.36111111111111111111111111111116"^ + 0.4254662447518595926183272929396 
- 31.041066530239171582155933528974e - 110.36756287612408836668594378878e^ 
611.019191920868811743653000126916^- 1734.0854215005636534710373849833e^ 

- 7316.7112322583252172396619274906e^ + O(e^) 



(6.28) 



Those values of Ii{d,n) are used to determine the initial condition for V5{d,t), Eg. s (|6. 1316. 1416.1^ . We 
then solve the differential equation (|6.5I) by expansions in series centered in the points i = 1, 1/2, 1/4, 
1/8, 1/16 and (as above, this subdivision is due to the presence of a singular point at t — —1/8). By 
inspecting the equation H6.5|l one gets that the behaviour at < = of V4{d,t) is 



V5id,t^ 0) 



.(1) 



{d)f 



-d+3 



„(2) 



id)t 



+ (d)t(-3<i+8)/2 ^ ^(4)(^)^(-5d+12)/2 



(6.29) 
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so that for c? — > 4 the integral (|6.4|l is surely convergent for n > 5; then we calculate the integral (|5.4|l 
for n — 5,6, 7, 8, 9 by integrating the series term by term. By using repeatedly "top-down" the recurrence 
relation Ht).16|l starting from n = 9, we obtain 15(0?, 6), I^id, 5) (used for cross-check), hid, 4), . . ., I^id, 1). 
The result, up to the coefficient of order 5 in (d — 4) included, is 



/5(d,l) =r(l + e) 



-0.25e~ 



0.875e" 



1.416666666666666666666666666667e" 



- 1.449652777777777777777777777778e"i - 14.055442461941065705599451765901 

- 9G.416062304531327135375791542063e - 1170.3684076603804614545918785105e2 
5299.6462727245240600241060624284e3- 37132. 219579420859394978093377604e'' 

- 144381.92488313453838475109116166e^ + O(e^) 



(6.30) 



Taking into account the normalization H6.2|l one finds that the first 8 terms of Eq. 1)6.30(1 agree with Eq. ((5.16|l . 

We want only to mention that we have also independently checked the numerical result ((6. 30(1 by 
calculating the master integral with all masses equal to one by difference equations, and then by using the 
value so obtained as initial condition for the integration of a differential equation in the photon mass A 
from A = 1 to A = 0. 



7 The a: ^ values. 

We evaluate in this section the x values of ^{d,x), Eo. 1(2. 11(1 . 
By combining the familiar formulae 



r(a + /3) } a;^-l(l-^)/3-l 



p - /)2 a)]"] {P + bf r(a)r(/3)7 [{I - xk)^ + x(l - x)k^ + ax + b{l - x)]°'+l^ ' 

and 



C-Ud) / 1 ^("-i) ^ (72) 

W / (2^)d_2 [(;_^)2 + c]„ 4 r(3-|)r(n) c"'- ' ^ ' 



one gets 



d'^l 1 1 1 r(a + /3-f) /• xi-P-\l 



-a-l 



{2^Y~^ [{k-iy + aY [P + bY 4r(a)r(/3)r(3-|) 7 [k^ + d{a,b,x)Y+P 
where 

ax + b{l-x) 

d{a,b,x) = — . 

x{l — X) 

We rewrite Eq. 1(2.11(1 as 

^ ' ' (2^)4('i-2) J [{p - qY + x2] j [q- k^Y J {k2 - kiY J [{ki - lY + 1](;2 + 1) ' ^ ■ ' 

and then use repeatedly Eg. ((7. 3(1 . using in the order the parameters y,yi,y2,z for integrating the loops 
I, ki, k2, q, obtaining 

m^) - r(i-2(d-4)) 



1024 {d ~ 3)(d - 4)(2d - 5)(2d - 7)r'i (l - \{d - 4)) 
1 1 1 

X 1 dy yi-^{l - y)i-^ j dy, yr\l - y,f-^ j dy^ yt^l y^Y 



X ^{d,x,y,yx,y2) , (7.5) 
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where 



nd,x,y,y.,y2) = Jdz . (7-6) 



and 

D{y,yi,y2) = -j-^ rr^ -r^ r • (7.7) 

y{l - y){l - yi){l - y2) 

The above formulae are vaHd for any x; form now on we take < a; <C 1. For definiteness, we take also d 

to be "just bigger" than 2 (i.e. d = 2 + rj, with < <C 1). The z-integral will be carried out first. To 
that aim, introduce an infinitesimal parameter Z, such that < a;'^ Z <C 1 (a possible choice might be 
Z = —x^ Inx, but the exact value of Z will be irrelevant), and split the integration interval as 



J dz = J dz + j dz 



z 

correspondingly, we write the z-integral as 

*(rf,a;,t/,yi,y2) = *i(d,a;,2/,j/i,j/2) + 1'2(d, a;, y, 2/1,2/2) , (7.8) 



z 



^i{d,x,y,yi,y2) = j dz 





1 



'^2{d,x,y,yi,y2) = j dz 



[{l-z)^x^+zD{y,y,,y2W-^'' 



[{l-z)^x^ + zD{y,y,,y2W-^'i " 
z 

In the second term we can neglect x"^ in the denominator obtaining simply 

1 „ 1 

'^^ T — m = D{y,yuy2f'^~^ / dzz^'^ 

[ zD{y,yi,y2) ]^ J 

z 

Diy,yi,y2f''-' , (7.9) 



z 

2 



d-2 

where we have neglected the contribution from the lower integration limit Z, which is Zi~^, as for d just 
bigger than 2 it vanishes with Z. 

The first term is slightly more delicate. To start with, asO<z<Z<^lwe can neglect z with respect 
to 1 in (1 - zf x^ 

15-2d ' 



for reasons which will be apparent in a moment, we rewrite it as 

a o 

22 ^ 



/■ z^-i"^ 

^i{d,x,y,yi,y2) ^ / dz -j-^- — — 

J [x^ + zD{y,yl,y2)]'■ 



'^l{d,x,y,yl,y2) = j dz ^2n5-2d 



P [Diy,y,,y2) + ^^ 
and integrate by parts the factor z^~'^; the result is 



^i{d,x,y,yi,y2) 



d-2] [x^+zD{y,yuy2W-^'' 

„2d-5 n f , z^-y 
= 2— — - x^ j dz 



z z 

2 



+ {2d - 5) a;M dz 




[a;2 + ^D(2/,yi,t/2)]6-2d " 



d-2 j [x^ + zD{y,yuy2)f-^'' ' 
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as the end-point contributions vanish for d just bigger than 2. We can now modify the integration interval 
for x, 2/2) into 

dz = J dz — J dz , 
z 

and write correspondingly 

'^i{d,x,y,yi,y2) = *ioo(d, a;, y, J/i, ^2) - yi, ^2) , (7.10) 

where, thanks to the previous integration by parts, the two resulting z-integrals (from to cx) for the first, 
from Z to cxD for the second) are now both convergent for d just bigger than 2. 
We start again from the second term. 



00 

J/2) = 2-^— y x dz 



[x^ + zD{y,yi,y2)]'^-^'^ ' 

z 

in the denominator we can neglect x'^, the resulting integral is trivial and the result can be written as 

^Md,x,y.y,,y2)^2-^ [i^(y, y„ ^2)]-- ^ (bj 
which vanishes with Z (recall x'^ Z), so that 

*iz(d,a;,7;,yi,y2) =0 . (7.11) 
In the first term of Ea. H7.10|l . 5'ioo(rf, 2;, y, j/i, 2/2), we substitute z = tx^ / D{y,yi,y2), obtaining 

vl/i^(d,x,y,2;i,y2) =2-— f [^(y, y2)]irf-4 / 



d-2 L"v«:^.,«.yj y - (H-t)6-2d 



By using the second representation of Euler's Beta function 

r(«)r(/3) 



1 00 

a-1 



du ^."-1 (1 - uf'' = / j^—y^ , (7.12) 



we finally obtain 



^ioo{d,x,y,yi,y2) = 2^^B(4-^d,2^^d] x^"-^'' [D{y,yi,y2)]i''-^ . (7.13) 



2d- 5 ^ / 3 , „ 1 
- — -S 4- -d,2 - -( 
2 V 2 2 

By collecting the results Ea.s H7.13l7.11l[7'.10l7.9l7.'5|l . one obtains the value of the function '^'{d, x, y, yi, ^2) 
to be substituted in Ea. H7.5|l : recalling Eg. 1)7.7(1 one sees that all the remaining integrations factorize and 
can be carried out in terms of Euler's Beta- functions Eg. ((7. 12(1 . 

It is clear that the final result for Eg. ((2. 11(1 for a; — > consists of two terms, the first constant (inde- 
pendent of x), the second proportional to x'^^'^. In the notation of Eg. ((3.5() we have 

4(1) ^ 3d -11 

° 8{d -2){d- 3){d - 4)3(2d - 5)(2rf - 7)(3d - 8)(3d - 10) 

r(i -{d- 4))r(i - 2{d - 4))r^ (i + ^{d - 4)) (i - l{d - 4)) 
^ r4(i-i(d-4))r(i-3(d-4)) 

,(2) ^ 2 (2d -7) 

° 3(d - 2)2(d - 3)(d - 4)4(3d - 8)(3d - 10) 

r (1 + i(d - 4)) r (1 - |(d - 4)) r2(i -{d- 4)) 



r2(i-i(d-4))r(i-2(d-4)) 



(7.14) 
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Figure 3: The 4-loop watermelon graph. 



Let us observe that the term Aq is the value of the vacuum graph in Fig. 3., in agreement with the result 
in Eq.(A.12) of [13] (up to a different normalization). 
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